In-Situ Dual-Port Polarization Contrast Imaging of Faraday 
Rotation in a High Optical Depth Ultracold 87 Rb Atomic 

Ensemble 

Franziska Kaminski * l , Nir S. Kampel 1 , Mads P. H. Steenstrup 1 , Axel Griesmaier 1 ' 2 , Eugene 

S. Polzik 1 and Jorg H. Muller n 

1 Niels-Bohr-Institute, Danish Quantum Optics Center QUANTOP, Copenhagen University, 

Blegdamsvej 17, 2100 Copenhagen, Denmark 
2 Physikalisches Institut, Universitat Stuttgart, Pfaffenwaldring 57, 70569 Stuttgart, Germany 

September 4, 2012 



Abstract 

We study the effects of high optical depth and density on the performance of a light-atom quantum interface. 
An in-situ imaging method, a dual-port polarization contrast technique, is presented. This technique is able to 
compensate for image distortions due to refraction. We propose our imaging method as a tool to characterize atomic 
ensembles for high capacity spatial multimode quantum memories. Ultracold dense inhomogeneous Rubidium 
samples are imaged and we find a resonant optical depth as high as 680 on the Dl line. The measurements are 
compared with light-atom interaction models based on Maxwell-Bloch equations. We find that an independent atom 
assumption is insufficient to explain our data and present corrections due to resonant dipole-dipole interactions. 

1 Introduction 

The storage and retrieval of single photons [ 1 1 and con- 
tinuous variable quantum states [2| in quantum memo- 
ries |3]|4) has become a major endeavor for the realiza- 
tion of quantum networks [ 5 1 . While single-qubit mem- 
ories are sufficient to establish a secure communication 
channel, being able to store more qubits increases the 
capacity of the channel. It has been shown that not only 
the fidelity of storage and read-out but also the multi- 
mode capacity of an ensemble scales favorably with the 
on-resonant optical depth (OD) |6|. 

Spatially resolved detection is a requirement for mul- 
timode memories [7] that store each polarization qubit 
in an independent spatial light mode. 

Polarization rotation, also called Faraday rotation, is 
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a well known effect [8]. It has been a means to gener- 
ate an entangled state of two atomic ensembles |9j, to 
measure the OD of atomic ensemble s |10||1 1) , spin dy- 
namics 1 12 13 1 and magnetic fields |14| and was used 
as a gigahertz-bandwidth probe (15). 

When applying the Faraday technique to dense, in- 
homogeneous, high OD samples new challenges arise. 
A large OD, together with the sample inhomogeneity, 
leads to stronger refraction or lensing, distorting im- 
ages. A small transverse size leads to diffraction, posing 
stringent constraints on the properties of the imaging 
system. A large density leads to effects beyond the in- 
dependent atom hypothesis commonly applied in quan- 
tum optics 1 16 17 1. 

We present an imaging method that reduces distor- 
tions due to refraction and we introduce a model that 
treats the influence of light assisted cold collisions on 
dispersive interactions. 
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Figure 1: Polarizability coefficients normalized to be 
unit free by 2|(7||ii||y')| 2 /^r / \. ff (0) describes the scalar 
refractive index, o- (1) is responsible for Faraday rotation 
and Q ,(2) is responsible for ellipticity changes and the 
corresponding Raman population transfers. The zero 
crossing of the coefficient is indicated as Ao and its 
maximum as A max . 

There are already several dispersive imaging tech- 
niques available. Phase-contrast imaging [18, 19) uses 
a spatially selective phase plate and is hardly polariza- 
tion sensitive. Therefore it mainly measures the scalar 
part of the polarizability, or scalar refractive index 
(Fig.[T|. Single-port polarization-contrast imaging [20] 
uses absorptive polarizers and detects the vector part, 
or refractive index differences but does not can- 
cel diffraction. Dark-ground imaging pT) uses a spa- 
tially selective block and is sensitive to both and 
a <l \ Another technique records directly the diffracted 
wave of an object and the image is reconstructed nu- 
merically 1 22 1 . 

We employ a polarization-contrast imaging method 
which uses a polarizing beam splitter (PBS) instead of 
an absorptive polarizer, enabling us to record the full 
light intensity (dual-port). This gives access to the 
scalar and vector components of the polarizability si- 
multaneously and distinguishably. It enables us to can- 
cel diffraction in the same way common-mode noise is 
canceled in a differential photodiode and leaves us sen- 
sitive to the spatial profile of the Faraday rotation signal. 
All the described methods are weakly sensitive to the 
tensor part of the polarizability, a (2) , leading to changes 
in the coherences and populations of the atoms during 
the interaction. This manifests in an excess ellipticity 




Figure 2: Optical setup showing the 45° PBS prepar- 
ing the probe light after the Fiber (F), the placement 
of lenses L\, L2,h and L^y, mirrors (M) and analyzing 
cube H/V PBS as well as the CCD camera. The B -field 
direction is mainly along z. The inset shows the 87 Rb 
D 1 line level scheme and the linear probe light in the z 
quantization axis. 

of the light. High sensitivity to the tensor components 
can be obtained with circular polarizers. 

A Raman/beam splitter-type memory J23) is based on 
the a <T> part of the polarizability making use of Raman 
population transfers and coherences. In our specific 
case of 87 Rb it can be realized using the nif — +1 states 
of the F=l manifold (Fig. |2j. To avoid unwanted dif- 
ferential phase imprints it is advantageous to minimize 
the a (1) part of the polarizability. On the Dl line a (1) is 
expected to vanish at a detuning Ao/(27r) = -204MHz 
(Fig. [I}. A local maximum of a m is expected for red 
detunings at A max l(2n) = -660MHz. 

Scalar diffraction is due to the a (0) part of the polar- 
izability that is dominant over the whole range of ex- 
plored detunings. An imaging method that can address 
these distortions is therefore desirable. 

Dual-port polarization contrast imaging, also referred 
to as Faraday imaging in this paper, can be more gen- 
erally employed for spatially resolved detection of the 
atomic magnetization and the study of quantum co- 
herence in degenerate gases and solid state systems 

2 Experimental Setup 

The ultracold thermal atomic ensemble is prepared in 
a Ioffe-Pritchard type magnetic trap with frequencies 
oj r = 2tt(1 15.4+0.5) s- 1 andw- = 2^(11.75+0.25) s" 1 . 
The atoms are spin polarized in the \F = \,nif = -1) 
state in the quantization axis defined by the local mag- 
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netic field. The main B-field component is oriented 
along the propagation direction of light (z). The gravi- 
tational sag of 18|xm leaves the B-field at an angle be- 
tween the z and gravity axis (-y) of about 15° in the 
center of the trap. There is a 5° variation within w r and 
a 0.25° variation within w z . 

The first element in the optical assembly, 
Fig. [2] is a polarizing beam splitter oriented 
at 45° to prepare a clean linear polarization. 
The probe beam enters the ensemble with an 
exp(-2) beam radius of 140[im with a flux of 
1100 photons /^s/(im 2 which corresponds to a sponta- 
neous emission probability of 0.03 in the pulse duration 
at a detuning AaI(2k) — -200MHz. The pulse duration 
is lO^is. To average over shot-to-shot atom number 
fluctuations we identically prepare each experimental 
run five times. 

After the atomic ensemble an achromatic lens dou- 
blet (Li) collimates the scattered light, which is then 
split into horizontal (H) and vertical (V) polarization 
components by another PBS and imaged with identical 
lenses Liy/H on t° separate areas of a CCD camera. We 
take images I v and I H with atoms present, images Iy l 
and F£ f without atoms present and bias images to cor- 
rect for any stray light and electronic offset. We balance 
the detection to split unscattered light to equal amounts 
into both output arms, Iy { = I^ f . Any small remaining 
imbalance is corrected for during post-processing. This 
allows us to detect the difference In - Iv and the sum 
Ih + Iv simultaneously with optimal sensitivity. Cen- 
tering the images Ih/v on top of each other is critical 
for the accurate determination of Faraday angles. We 
therefore apply a fitting algorithm on images at large 
detuning and the determined position of the images is 
then used for the whole data set. We achieve sub-pixel 
resolution for the centering. 

The Faraday angle Of is the rotation angle of linearly 
polarized light. We deduce the angle on each pixel by: 

a 1 . ( Ih-Iv \ m 
6 F = - arcsin . (1) 

^ \Ih + Iv! 

This will give an accurate polarization rotation angle, 
if we can neglect the presence of any circular polariza- 
tions after the interaction. The sensitivity of our method 
is illustrated by our ability to detect a small rotation of 
0.038°, produced by the cell windows that are subject 
to large magnetic fields due to the magnetic trap, for 
which we correct the data. 

We adjust L\ to image the end plane of the ensem- 



ble. The diffraction limited imaging resolution on the 
object plane is 3.6|im. The magnifications with 95% 
confidence bound in both arms of our imaging system 
are 12.83 ± 0.08 and 12.87 ± 0.08, i.e. identical within 
measurement accuracy. Images of clouds falling freely 
under the influence of gravity are used to determine the 
magnification. The acceptance full opening angle of the 
system is 15.2°. 

For quantitative imaging the opening angle of the 
imaging system needs to be sufficiently large compared 
to diffraction angles of the ensemble. We estimate 
the full geometric diffraction angle of the ensemble is 
arc ~ A I Ad < 1°, where A is the wavelength and d the 
radial extent of the ensemble. The angle due to refrac- 
tion or lensing can be approximated by comparing the 
phase shift <p in the center to the one at the edge of the 
sample a L = 2((f>(r = 0) - <p(r = d))A/(nd) (21) and 
reaches at A A /(2n) = -200MHz a L (A ) * 0.3°. Both 
are sufficiently small to neglect light loss at the aperture 
of the imaging system. We compared these estimates to 
more elaborate diffraction simulations [6 25| including 
a model of the full imaging system. These confirmed 
the conclusions of the simple diffraction estimates. 

As an independent sample characterization we per- 
form standard absorption imaging on the D2 line 
after a time of flight (TOF) of 45 ms. Us- 
ing the magnification 1.577 ± 0.002 and indepen- 
dently measured trap frequencies we can spec- 
ify the waists of the in-trap density distribution 
p(r,z) = poexp(-r 2 /(2w 2 ) - z 2 /(2w 2 )) as w r = 7.7um 
and w z - 70.5 um, leading to a Fresnel number of 
F = w 2 r l(2w z A) = 0.5. From these parameters we 
determine the temperature k B Tj = Mw 2 to 1 and T = 
(TyT z ) ' 3 = 300nK, where M is the atomic mass. 

By using the D2 line scattering cross section for the 
cycling transition ctdi - 3/l 2 32 /27r we determine a peak 
density of p? s = 1.2 • 10 19 m -3 and an atom number of 
N^ s = 8.1 • 10 5 as an average over all data points. 

Both absorption imaging and Faraday measurements 
allow to deduce the number of atoms using models for 
the optical cross sections. Comparing the deduced num- 
bers and taking into account the thermodynamic prop- 
erties of ultracold gases allows to identify systematic 
errors in either method. While absorption imaging is 
a standard method, it is well known that it is difficult 
to estimate the precise effective scattering cross section 
due to uncertainties in magnetic field alignment, light 
polarization quality and repump efficiency 1 26 1 . Since 

In later control experiments we found a systematic undercount 
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the cycling transition allows for the maximal cross sec- 
tion ctd2 our measured A^ t bs and py bs are hard lower 
bounds. We estimate a hard upper bound by noting that 
we do not observe condensed atoms on the absorption 
images and hence T/T c > 1. The condensation tem- 
perature of an ideal gas is ksT c - hu>(N a ,/£(3)) 1 ^, with 
a> = (a>^(L> z ) 1 ^ and £(3) is the Riemann zeta-function. 
We correct T c for the effects of finite size and mean field 
interactions. Both effects reduce T c for our parameters 
by maximal 1.2% and 5.5% respectively |26]27). Cor- 
rections to the determined ensemble temperature T arise 
from the use of a Gaussian instead of a Bose distribution 
in the fit model. We estimate this systematic correction 
by fitting Gaussian profiles to analytical Bose enhanced 
densities and find a systematic underestimation of the 
temperature by about 10%. 

Taking all these corrections into account we reach 
T/T c = 1 for several single-run data points at an atom 
number scaling factor / = N a i/N^ s of maximal f max 
2.7 . defining a hard upper bound for the real atom num- 
ber N ar . We will show below that this upper bound is 
still too low to allow the Faraday rotation data to be fit- 
ted with an independent atom model and we conclude 
that lineshape corrections due to light assisted collision 
become relevant. 



3 Experimental Results 

Figure [3] shows the detuning dependence of the ob- 
served peak Faraday angles (black squares) together 
with a model including light assisted cold collisions 
(red line) and a coupled Maxwell-Bloch model that as- 
sumes independent atoms (grey area). Both models are 
discussed in Sec. [4] Both use input parameters de- 
duced from absorption imaging, i.e. the sample radii 
w r and w z and the atom number fN^ s , averaged over 
all detunings. We infer an optimal atom number scaling 
/, by matching the light assisted cold collision model 
to the experimental data and obtain f opl = 2.13 (0). 
The grey shaded area indicates the atom number scal- 
ing range 1 < / < 2.7, defined in Sec. [2] Experimental 
peak angles are determined by averaging over 3x3 pix- 
els around the determined center positions of the den- 
sity distribution and the error bars are the standard devi- 
ation of 5 experimental runs. The small structure on the 
red line is due to vibrational molecular resonances, dis- 
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of atoms by a factor of two in absorption imaging due to insufficient 
repumping. 

2 This corresponds to a temperature ratio of T/T c = 1.27. 
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Figure 3: Detuning dependence of peak Faraday angle. 
Experimental data (black squares) is best reproduced 
by a model including light assisted cold collisions (red 
solid line). The grey area shows the prediction of a 
coupled Maxwell-Bloch model assuming independent 
atoms for the permissible range of the atom number 
scaling factor / (see text). The inset shows the relative 
atomic density variation with detuning. 



cussed further below. The figure inset shows the vari- 
ation of densities with detuning normalized to the av- 
eraged density entering the models. The discrepancy 
between the data point at A A /(2n) = -1340MHz and 
the models might be explained by the low density. 

Figure |4] shows the spatially resolved Faraday angle 
as deduced from the camera images for the detunings 
A A /{2n) = -{1000, 230, 200} MHz averaged over five 
realizations together with the two model predictions for 
input atom number f ptN& = 1.73 • 10 6 . The good re- 
producibility of sample preparation is evidenced by the 
small standard deviation encoded in the colored areas 
around the averaged profiles. The spatial shape of the 
Faraday angle profile at A A /(2n) - -1000MHz (red) 
fits the expected shape from the light assisted cold col- 
lision model when the finite imaging resolution is taken 
into account. The experimental profile at A A l{2n) = 
-200MHz deviates significantly from the shape of the 
density distribution. We observe minimal Faraday rota- 
tion at A A /(2n) = -230MHz, shifted by about 30MHz 
from the expected detuning Ao. This shift is outside 
possible systematic errors in the frequency scale. 

Figure [5] shows a three pixel averaged cut through 
the transmission T = (/# + 7y)/(/J| f + /^ f ), where we 
normalize with the reference images. This allows us to 



4 




-40 -20 

In Trap Y Coordinate [|im] 

Figure 4: Cut through Faraday angle images for var- 
ious detunings At. Colored areas represent the stan- 
dard deviation of 5 experimental realizations. Cou- 
pled Maxwell-Bloch simulations assuming independent 
atoms (squares, circles) underestimate the angle by a 
factor 1.46, compared to a model including light as- 
sisted cold collisions (dashed, dot-dashed). 
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Figure 5: Spatially resolved transmission: cut through 
T = (Ju+IvMir^ +Iy f ) for various detunings compared 
to the column density approximated transmission pro- 
file of the ensemble at A A /(2n) = -200MHz (dashed). 
Diffraction dominates over absorption for all detunings. 



visualize the effect of intensity redistribution across the 
image due to refraction and diffraction. To indicate the 
expected photon loss we plot an estimated transmission 
profile for A^/(27r) = -200MHz using a naive column 
density model neglecting diffraction effects. The ex- 
pected photon loss is hardly distinguishable from the 
detection noise. From the spatial transmission curves it 
is apparent that data is dominated by refraction rather 
than absorption for all detuning values shown. We em- 
phasize that due to the dual-port detection the distortion 
effects of diffraction are largely canceled in Faraday an- 
gle profiles. 

This compensation of refraction effects is illustrated 
in Fig. |6j where raw images In, Iv, F^f an d are 
shown together with the 2-D reconstruction of Faraday 
rotation angles measured at a detuning of A / i/(2^) = 
-400MHz. 



4 Models 

Our first model accounts for the collective response of 
all atoms, while treating each atom as an independent 
scatterer. The model is based on coupled Maxwell- 
Bloch (MB) equations |23]|28]|29) with excited states 
eliminated adiabatically, using continuous variables and 
including absorption. The probe light is propagated 
spatially through the ensemble while simultaneously 
evolving the ground state F= 1 manifold populations and 
coherences in time. To incorporate the spatial inhomo- 
geneity of sample density, initial atomic state and mag- 
netic field, we extend the 1 + 1 dimensional geometry to 



3+1 dimensions, following |30|. In the model light is 
assumed to propagate along straight lines and atomic 
motion is neglected. 

The model Hamiltonian <H = tf ] + e HP r +'H (2 ] +"K B 

ml ml mt D 

contains the atom-light interaction decomposed into its 
irreducible tensor components < H^ (App. |Xj) and the 
effect of the external magnetic field "Hb = • F. Here 
Q.(r) is the the vector of Larmor frequencies and F is the 
total atomic angular momentum vector. To compare the 
simulation results to our image data we plot the Faraday 
angle, time averaged over the probe pulse duration, at 
the output end of the atomic sample (Fig.[3]and|4]i. 

We can use the full MB model to quantify the com- 
bined effects of tensor polarizability and B-field inho- 
mogeneity by comparison to a much simpler, idealized 
Faraday model. 

The idealized Faraday model starts from an effective 
two-level system description extended to our multilevel 
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Figure 6: Raw images (left panel) and deduced Fara- 
day image (right panel) for a detuning of Aa/(2tt) = 
-400MHz. V/H refers to images with atoms present, 
while V/H Ref to images without atoms. The color 
scale indicates the Faraday angle in degree. The scale 
bar indicates dimensions on the camera. Diffraction 
rings on raw images V/H disappear on the Faraday im- 
age. 

situation and neglects atomic evolution due to rank 2 
tensor components and the magnetic field. The Faraday 
angle in this model is half of the accumulated phase dif- 
ference between the left and right circularly polarized 
light. The phase shift (p for a single polarization com- 
ponent is given by </> = 2nA~ 1 J ndz, where n is the re- 
fractive index and A is the wavelength. Summing over 
all excited states (z) we obtain a Faraday angle 



0£(Aa) 



^p'^M^jp^ (2) 



where is +1 for right-handed and -1 for left- 
handed circular polarization, p(r) is the atomic density, 
Fa is the full atomic line width and A® is the detun- 
ing to the excited state i. The scattering cross sec- 
tion is cr®(A®) = of /(I + 4(A%/T A ) 2 ) with cr« = 
■ 3A 2 /(2tt) ■ (27' + l)/(2J + 1) the on-resonant scat- 
tering cross section. Here J and J' denote the total 
electronic angular momentum of the ground and ex- 
cited state respectively. The parameter is defined via 
the transition dipole moment d, = gi(J\\d\\J'), such that 
in a two level system £ = 1. Comparing this simpler 
model to the full Hamiltonian dynamics we find that 
the effects of inhomogeneous magnetic field and tensor 
dynamics lead to a reduction of the Faraday angle by 
a constant factor /3g = 0.86 in the range of detunings 
A A /(2n) = -1340MHz to -400MHz. 

Our second model addresses the effect of light as- 
sisted cold collisions. At high atomic densities atoms 



can no longer be treated as independent scatterers. Elec- 
tronic energy levels for close pairs of atoms split and 
shift. The light scattering properties of a pair are mod- 
ified compared to isolated atoms. This effect of the 
dipole-dipole interaction can be described by estab- 
lished methods from molecular physics 1 3 1 32] . We 
consider repulsive and attractive molecular potentials 
for ground-excited state Rb* 2 atom pairs, neglecting hy- 
perfine recoupling p3), and calculate allowed energy 
levels. For the attractive molecular potentials the posi- 
tion of photoassociation resonances are calculated using 
the LeRoy-Bernstein formula |34|. On repulsive poten- 
tials atom pairs can be excited to a continuum of states. 
We are interested in the dispersive effects of all these 
shifted optical resonances. Writing the additional Fara- 
day rotation analogous to the simplified Faraday model 
we obtair0 



1 a® r°° 

er(AA) = iJ> r7 J ^fMA A ) P ^dz. (3) 



Here Ay 1 ' is the detuning to the molecular resonance 
v for a specific excited state i, F v is the linewidth 
of that resonance. The pair absorption cross section 
<>pa 4a) depends on the density of pairs at a given 
detuning A^. Via the definition cr scat = <£> SC atjB, where 
I lnc is the incident intensity and hco is the energy of a 
photon, we can relate the pair absorption cross section 
to the scattered photon flux <& scat . The scattered photon 
flux is identical to the atom pair loss rates given by Juli- 
enne et al. [35 36], since exactly one photon is scattered 
per atom pair loss event (see App. |E) . 

The pair absorption cross section crpa is linear in 
density (App. [B]i and therefore the light assisted cold 
collision Faraday angle is quadratic in density. Thus 
the spatial profile of the pair contribution to the Fara- 
day angle has a reduced Gaussian width W r a = w r /V2. 

We calculate the total Faraday rotation angle by using 
equations [2] and [3] To correct for magnetic field inho- 
mogeneities and tensor evolution we multiply the result 
by the above defined /3b, such that 



(4) 



The total optical depth OD(A A ), defined via the inten- 
sity attenuation I/Io = exp(-OD{A A )) is given by the 
product of the atomic density and the scattering cross 



F. Kaminski, N. Kampel, A. Griesmaier, E. Polzik, and Jorg H. 
Miiller, to be published 
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Figure 7: Calculated absorptivity and Faraday angle en- 
hancement due to interactions between atoms. Plot- 
ted is (OD* + ODP") IOD A and + 0^) /6* for an 
atomic density of foptP^ = 2.6 • 10 I9 rrT 3 . The inset 
shows the detuning and hyperfine structure corrected in- 
dependent atom optical depth. 



section integrated along the propagation direction of the 
light: 



= OD A (f ± ,A A ) + OD^ ±< A A ). 



(5) 



From the determined atom number f op tN^ s and en- 
semble size w z we calculate the independent atom on- 
resonant OD for the |1,-1) to |2, — 2) Dl transition 
with i* = 1/2. We find a peak optical depth OD A = 

g(3A 2 m /2n)p A& z = 680. 

We present the predicted enhancement due to 
light assisted cold collisions of absorptivity (OD^ + 
ODP")/OD A and Faraday angle (Q A + QP a )jQ A for the de- 
termined atom number as a function of probe detuning 
Aa in Fig. [7] The inset gives the independent atom OD, 
taking into account all relevant levels for Faraday rota- 
tion (see Fig.|2]). The enhancement (OD A + OD pa )/OD A 
of the absorptivity due to atomic interactions is about 
25% in between resonances and even up to a factor of 
10 close to vibrational resonances in the plotted detun- 
ing range. For the Faraday angle we find an almost con- 
stant enhancement of 1.46 at all detunings away from 
Ao. The local enhancement due to single vibrational 
resonances is small. 



5 Discussion and Conclusion 

The data presented in Fig. |4j [5] and [6] show that the 
influence of diffraction on Faraday images is reduced. 
Common-mode diffraction on both images, ly and 
is canceled when calculating the difference In - /y. The 
common-mode diffraction stems from the scalar polar- 
izability a' ', which is the largest contribution to the 
polarizability as shown in Fig. [T] This compensation 
is, however, not perfect and we discuss in the following 
the effect of uncompensated diffraction and refraction 
on the detected Faraday angle. In a geometric optics 
picture, the trajectories of light rays are curved due to 
the inhomogeneous density of the sample. This leads 
for strong refraction to a breakdown of the column den- 
sity approximation, which implicitly assumes a straight 
line ray path. Eventually this leads to noticeable differ- 
ences between the column density and the spatial profile 
of the Faraday angle. For red detunings the extended 
atomic cloud acts as a thick collimating lens, such that 
ray trajectories are bent towards the center of the cloud, 
which leads to a reduced Faraday angle in the center. 

Differential diffraction and lensing is associated with 
the a m part of the polarizability. It leads to a mismatch 
in wavefront of the left- and right-handed circular polar- 
ization modes at the exit plane of the atomic ensemble. 
This introduces locally ellipticity to the initially pure 
linear polarization. Since detection in the H/V basis is 
insensitive to circular polarization this lowers the Fara- 
day angle by a second order correction. 

While residual diffraction and refraction reduces the 
Faraday angle, Fig. [3] shows measured peak Fara- 
day angles significantly above the prediction from the 
Maxwell-Bloch simulation, which assumes the atoms to 
be independent scatterers and does not include diffrac- 
tion. We match our light assisted collision model to 
the data by choosing f op , = 2.13 < f ma x, scaling 
our inferred atom number from absorption imaging to 
f opt N% s = 1.73 10 6 . Comparing the cold collision 
model for this input atom number to the corresponding 
Maxwell-Bloch simulation we find an increase of the 
Faraday angle of 1 .46 as indicated in Fig. [4] Trying to fit 
the data directly with the Maxwell-Bloch model would 
require an atom number scaling of \A6f opt = 3.1. This 
lies 15% above the conservative upper bound of f max = 
2.7, discussed at the end of Sec.|2]and strongly suggests 
that the independent scatterer assumption breaks down. 

Our strategy to correct the optical response of the gas 
by considering atom-atom interaction on molecular po- 
tentials and the corresponding redistribution of oscil- 
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lator strength in frequency space can be contrasted to 
other approaches to describe the optical properties of a 
dense gas. Instead of calculating the collective response 
by a systematic expansion in atomic density ]37| or by 
a configuration average over many randomly placed in- 
teracting point dipoles p8[ , we focus on the contribu- 
tion of close pairs for which big resonance shifts oc- 
cur which lead to important modifications in the wing 
of the atomic line. The number of close pairs is deter- 
mined using the quantum mechanical scattering wave 
function for atom pairs interacting on the ground state 
molecular potential, hence particle correlations are ac- 
counted for. Close to the unperturbed atomic resonance 
we expect our approach to fail, since for the correspond- 
ing large internuclear distances pairs can no longer be 
considered isolated and the collective response of ever 
bigger clusters of atoms should be calculated instead. 
Our model predicts a surprisingly large modification of 
the Faraday rotation angle even for the modest particle 
density p2* s / opt A 3 = 0.05 used in the experiment. The 
model neglects hyperfine recoupling on molecular po- 
tentials and does not explain the observed shift in the 
position of Aq. Interestingly, in a recent experiment, 
which employs resonant absorption imaging as a detec- 
tion method for high density 2-D quantum gases fT6| , 
a decrease of absorptivity with increasing density has 
been observed. This is consistent with our simple pic- 
ture of redistribution of oscillator strength from the line 
center into the wings due to the resonant dipole-dipole 
interaction. 



We turn to the suitability of our atomic samples to 
multimode quantum memories. With the favorable op- 
tical depth and Fresnel number a mode capacity in the 
hundreds is predicted in forward read-out p9| . In ex- 
perimental implementations this number will likely be 
limited by the finite resolution of the imaging system. 



While the increased Faraday angle signals a higher 
coupling between atoms and light it remains to be seen 
in future experiments to what extend decoherence is in- 
creased by the resonant dipole-dipole interaction. For 
this, the presented weakly destructive dual-port detec- 
tion method will be an invaluable tool since details of 
the radial spin density distribution can be examined re- 
peatedly in face of strong refraction and diffraction ef- 
fects. 
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A Light-Matter Interface Simula- 
tion 

The interaction Hamiltonian in continuous variables for 
light propagation along the quantization axis is: 

( HZ ^H u a (0 \A A )Sot (6) 

<H% =H a m (A A )S 3 P z (7) 

<H% =H a (2) (A A )x 

[S , (P 2 X - Pfj + S 2 (P x F y + P y P x ) - 
§ (P 2 z -lF(F + V)l)], (8) 

where we defined 

and omitted the space and time (z,t)-dependence of 
all operators. A is the interaction area, d A is the transi- 
tion dipole moment and eo is the vacuum permittivity. 

The explicit commutation relations for continuous 
variables are: 

[Si(z,t),Sj(z',t')] = i^eijkSkiz^Wz-z'W-t') 

k 

(10) 

[Mz, t), fj(z', 0] = \ Yj e '*/*fc - z')S(t - O 

" k 

(ID 

and the Stokes operators with units [S] = [N p h/cT] = 
1 /meter, number of photons N p h over speed of light c 
and interaction time T are defined as, again omitting 
(z,t): 
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I-45) 



(12) 
(13) 
(14) 
(15) 



We solve the equations of motion for the F=l man- 
ifold with all three ground state levels. We choose to 
propagate the density matrix & and the second order co- 
herence matrix for the light p - (S 1 + §2 + S{)l%: 



dp 

a& 

~dt 



=i[p,n n ,+'H B ] 

=i[&,<H ira + < H B \. 



(16) 
(17) 



The Stokes operators are related to the experimen- 
tally taken images by S 1 oc I H — I v and So oc I H + I v . 
The Faraday angle Of is the rotation angle of linearly 
polarized light and is thus half the angle of the Stokes 
vector rotation in the S1/S2 plane of the Poincare sphere 
and reads in Stokes parameters 



n 1 • (Si 

9f = - arcsin — 
2 \S 



(18) 



B Light Assisted Cold Collision 
Parameters 

The partial atom loss rate due to photoassociation to the 
level (z, v) is given by [35 36 1 



?& v V,Aa) 



2p(r)K v g 



v v r v 



A« 2 + (r v /2) 2, 



(19) 



where T v is the full linewidth of the molecular res- 
onance v and Ay is the detuning relative to this reso- 
nance. The parameter v v characterizes the slope of the 
excited state molecular potential in frequency units and 
K v is the rate coefficient. 

As explained in Sec. [4] the partial pair absorption 
cross section is then given by 



y^ v) (r,A A ) hcj 



(20) 
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